Circuit QED bright source for chiral entangled light based on dissipation 
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Based on a circuit QED qubit-cavity array a source of two-mode entangled microwave radiation 
is designed. Our scheme is rooted in the combination of external driving, collective phenomena 
and dissipation. On top of that the reflexion symmetry is broken via external driving permitting 
the appearance of chiral emission. Our findings go beyond the applications and are relevant for 
fundamental physics, since we show how to implement quantum lattice models exhibiting criticality 
driven by dissipation. 



The controlled generation of quantum states of light is 
crucial for applications both in quantum information pro- 
cessing [1] and precision measurements [2] . While Gaus- 
sian entangled states of light can be generated through 
non-linear optical interactions [3], the field of supercon- 
ducting quantum circuits and circuit- QED provides us 
with more efficient tools to explore such applications 
[4]. In this respect, we recall the preparation of single 
mode squeezing [5], two-mode squeezing [6, 7], broad- 
band squeezed light [8] single-mode nonclassical states of 
light [9] and entanglement [10], as well as the extremely 
precise tomography both in cavities [9] and in open trans- 
mission lines [11]. 

Two distinctive features of circuit-QED (cQED) are 
the precise positioning of quantum emitters and their lo- 
cal control through driving fields or, more recently, dis- 
sipation [12]. These ingredients enable a paradigm shift 
in the control of light-matter interaction. First, a con- 
trolled separation between quantum emitters allows us 
to tailor collective effects like directionality and interfer- 
ence, with an accuracy that can hardly be matched by 
atomic systems. Second, these collective effects can now 
be embodied with dissipative dynamics [13], exploring 
novel dissipative quantum many-body phenomena [14- 
16], such as non-equilibrium phase transitions and dissi- 
pative criticality classes. In this work we combine both 
ideas, presenting a source of entangled light where the 
directionality and frequency of the entangled modes is 
tuned by means of collective phenomena, and where the 
source of entanglement is a dissipative many-body sys- 
tem that approaches criticality. 

The device that we present is a quantum metamate- 
rial that embeds a regular array of dissipative quantum 
emitters, merging a variety of physical effects into a sin- 
gle cQED device: (a) Non-perturbative sideband transi- 
tions in the qubit-cavity couplings, engineered by suit- 
able qubit drivings, (b) An incoherent qubit re-pumping 
mechanism based on a combination of local reservoirs 
with the qubit driving, (c) Directional emission of light 
in the cavity array supported both by a collective inter- 



ference and a the band engineering in this metamaterial. 
Through a detailed theoretical study of this setup we 
prove that our device has the following properties: (i) It 
acts as a source of entangled Hght, created by a dissipa- 
tive process that arise from a combination of precisely 
tuned drivings with fast qubit decays, (ii) The degrees 
of entanglement and squeezing, as well as the frequency 
of the entangled modes, can be tuned by means of the 
periodic drivings, (iii) It implements a chiral dissipa- 
tive process with broken reflexion symmetry, induced by 
the phase pattern of the driving. In practice, this ef- 
fect allows us to control the direction and momentum 
of the emitted entangled modes, (iv) The setup can be 
used to study dissipative many-body collective effects, 
both because of its non-trivial collective steady-state, but 
also because the system undergoes a dissipative quantum 
phase transition for suitable parameters. 

The setup.- Our dissipative source of chiral quantum 
light is schematically drawn in Fig. la. It consists on a 
ID arrangement of coupled cavity-qubit units where the 
qubits gaps are periodically and uniformly driven, and 
where each unit is coupled to an independent thermal 
bath to which they can dissipate. We describe the coher- 
ent part of the cavity-qubit units using the Hamiltonian 

^CQ = (^r^a, + ^a] + ga^Xj^ + J ^ XjX, (1) 

j {id) 

Here cr| labels a qubit state with frequency e, Oj rep- 
resents the creation of a photon with frequency cj^, 
and Xj = [aj + a^). The Hamiltonian includes both 
cavity-qubit, ^, and inter-cavity coupling, J. For mod- 
erate couplings and no resonant driving, the cavity- 
qubit and cavity-cavity interactions can be modeled 
in the Rotating Wave Approximation (RWA), form- 
ing the so called Jaynes-Cummings lattice (JCL), with 
^(^/% + H.c.) + Y.^..^ Ji^cij + H.C.). Analogously 
to the Bose-Hubbard model, it supports an insulator- 
superfluid quantum phase transition [17-22], which orig- 
inates from the competition between the repulsive non- 
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FIG. 1. (a) Setup of the system studied. Each supercon- 
ducting resonator is coupled to a qubit. The resonators are 
coupled each other, e.g. capacitely. We also drawn arrows 
representing the classical driving and the losses, (b) Sketch 
for the quantum light source. The coupled cavity can be em- 
bedded inside an open transmission line. The entanglement 
generated in the array is mapped onto the line. 



linearity induced by the two-level system and the hopping 
term J. A promising field to implement such many-body 
systems is cQED, as already chips with dozens of coupled 
identical resonators have been reported [23, 24]. 

In contrast with ordinary JCL's, for us the qubits do 
not act as a source of nonlinearity, but rather mediate an 
engineered dissipation. Furthermore, our model strongly 
relies on the use and control of external qubit drivings. 
All together the dynamics is given by a Linblad-type mas- 
ter equation [25] that reads 



dtQ = - i[HcQ ^Hacq] 



{cr+cr^- , g} 



(2) 



Through the paper, we neglect radiative losses, which 
can be straightforwardly included in our formalism and 
do not modify qualitatively the discussion below. The 
first important ingredient in this model is a uniform two- 
tone driving of all qubit gaps 



HaC = ^ COs{nrit + ^n,j)crj 

n=l,2 j 



(3) 



with tunable phases, 0i,2j- As we will see this driving 
plays a major role in the collective effects. The second 
important ingredient is the cavity and qubit dissipations, 
K, and r. In this letter we work in the limit of strong qubit 
decay, F ^, in which the quantum emitter customizes 
the effective environment of the photons. 

Dissipation engineering.- Dissipation is normally re- 
garded as the worst enemy for preserving quantum coher- 
ence in general and entanglement in particular. However, 



an appropriately engineered dissipation is a very efficient 
and robust way to drive a system to the desired quantum 
state. This powerful idea has been considered not only 
on single-particle models such as laser cooling [26, 27], 
but also more recently in the engineering of strongly cor- 
related phases [28] and distant entangled states [29, 30]. 

The principle underlying all these examples is a tai- 
lored system-bath interaction, which allows us to engi- 
neer the dissipator C in the master equation governing 
the irreversible dynamics, dtp = Cp. These dissipators 
drive the system of interest to its stationary state, >C^* = 
0. In most relevant cases, Davies' theory assures conver- 
gence to g* for any initial condition [31], thus avoiding 
the need to initialize quantum states. All this together 
makes dissipation engineering a promising new paradigm 
for quantum information processing [14, 15, 32]. 

Our source of entangled light is also based on engi- 
neered dissipation. In absence of driving, the setup from 
Fig. 1 thermalizes via Eq. (2) to a photon vacuum. How- 
ever, as soon as we modulate the qubit energy levels, the 
effective dissipation changes its asymptotic state, p*, to 
a product of two-mode entangled states. To make our 
arguments clear, we will first sketch the main ideas using 
a single qubit-cavity system and then introduce the pho- 
ton hopping, J, leaving all details for the supplementary 
material [33]. 

Following Refs. [26, 27, 34], we prove in three steps 
that a bad qubit can cool a cavity to a squeezed vac- 
uum. First, we adopt an interaction picture with re- 
spect to the qubit, resonator and driving, i^cQi = 
g{ae-'^^^ + h.c.)(cr+ e^-^^^) + h.c), where f {t) = et ^ 
^a=i 2i^(^/^(^) ^^^{^at-^4'a) iucludcs the periodic driv- 
ing (3). Next, under the assumption of weak driving, we 
use the Jacobi- Anger expansion [35], retaining terms up 
to order (Xa/fta)- The result in the RWA is 
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h.c.) 



(4) 



for driving frequencies 1^1,2 = A + e±cj^. Note that the 
qubit is now coupled to the squeezed resonator modes 



6 = a + rje'^^a^ 
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(5) 



Here and in the following we consider for simplic- 
ity unnormalized squeezed modes: [b^b^] ^ 1. The 
squeezing and the effective coupling strength, g = g x 
are determined by the external driving, gi = 
Ji(2AiM)Jo(2A2M) and g2 = Jo(2AiM) Ji(2A2M), 
through Bessel functions. The final step consists in 
an adiabatic elimination following the hierarchy of time 
scales e F uOr^g^/S.. The first inequality vali- 
dates Eq. (2) since it allows to treat the qubit dissipa- 
tion in a weak-coupling or Lindblad type master equa- 
tion. The second inequality justifies not only the RWA 
in the previous steps, but it also indicates that the 
qubit relaxation, F, is faster than its interaction with 
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the cavity, and it may be thus adiabaticahy ehminated. 
The result of these manipulations is a master equation, 
dtQ = ^^{2bgb^ — {b^b^g})^ that cools the resonator to 
the vacuum of the squeezed mode, b. 

What happens when the on-site squeezing of each res- 
onator competes with a coupling between resonators? 
Since the dynamics of both processes happens on dif- 
ferent operator basis, one would a expect a competi- 
tion between these phenomena and even some phase 
transition. For answering this problem we now move 
on to the coupled Hamiltonian (1). In the simplest 
case of translational invariance and periodic boundary 
conditions, we can introduce momentum space modes 
ak = A^"^/^ e"^^%^- {k e 2ti/N x Z), and show that 
the problem is similar to the single resonator case. More 
precisely, following the same definitions and approxima- 
tions as before, we obtain the effective master equation in 
the appropiate interaction picture (see Supp. Mat. [33]) 



2g 

dtQ = ^ -iujk[a\ak, q] + —{2bkQb\ - {b\bk, q}). (6) 



The chain of qubits is now cooling the resonators to 
the vacuum of two-mode-squeezed operators b^ = aj^ -\- 
V(^^-k-\-q ^^^^ dispersion relation (J <C Ur) Uk = A -\- 
2Jcos(/c). Let us remark how the external driving fully 
determines the properties of the asymptotic state. In 
particular, while the phase of the driving, <7 = (/>i — (/>2, 
selects the pairing between modes in we will show 
that the choice of frequencies, A = {fti -\- — 2e)/2, 
customizes the band structure and the amount of entan- 
glement. These are the main practical results in this 
manuscript. 

Entanglement in the stationary solution. — After the 
adiabatic elimination, the effective master equation (6) 
can be written as a direct sum of quadratic dissipators 
acting on the Fock spaces of the operators {bk^b-k+q). 
Consequently, the asymptotic state of the master equa- 
tion will be a product of Gaussian states in each of 
these Hilbert spaces, g* = ^p^-k+q- Each of the fi- 
nal density matrices -k+q ^^^^ f^l^Y characterized by 
the first and second moments of the operators R^^^ = 
{Qk,Pk,Q-k+q,P-k+q}, with Qk = l/V2{al + a/e) and 
Pk = —i/V2{al. — ak)' In particular, the first moments 
are all zero, while the second moments are conveniently 
grouped in the the two mode covariance matrix 

% = \{nl''R'Jl' + = ■ (7) 



This features a diagonal matrix a = {1- 
and two nonzero off-diagonal blocks 



-r?2)/2(l-r?2)i2, 
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-k+q 



(8) 




FIG. 2. (a)Logarithmic Negativity En from Eq. (10). (b) 
Cuts of (a) at the values of ^ = 0,0.2 and 0.4 (lower), as 
arrow-marked in (a). For values of 8 above 0.5 becomes 
negligible. 



Note that the most relevant parameter is 



^k-k+q = r(cjfe 

= F(A- 



hcj_fe+g)/4^i 

• 2J cos{q/2) cos(/c 



(9) 



because it determines the degree of entanglement be- 
tween the pairs of modes, {k^ —k -\- q). We quantify this 
through the Logarithmic Negativity [33], which is 



En = - log2 



1 + r 



2r] 



^k-k+q 



+ (1_^2)2 



(10) 



whenever Ek-k+q < (1 ~ rfY^'^ /r] and zero elsewhere. 

Fig. 2 shows a contour plot of En in terms of both 
£k,-k+q and 7^, together with some sections for fixed 
8. Let us remark how entanglement significantly grows 
when 8k^-k+q approaches and when r] approaches 1. 
A qualitative explanation follows from the master equa- 
tion (6): when £k-k+q approaches zero, it means that 
uok + ^-k+q = 0- Consequently, the terms that generate 
entanglement, a\po}_^j^^-{a\a^_^j^^, g) and 2akga-k+q- 
{aka-k+q-, q}i do not oscillate and are not suppressed. 
Moreover, since the strength of those terms grows as 
rf' and we know that 77 < 1 by definition, the optimal 
amount of entanglement is found when r] approaches 1. 
We focus on the case < 1, which corresponds to cooling 
to a squeezed vacuum, since the case 77 > 1 does not have 
a well defined steady-state, thus signaling an instability 
in the system. 

We finish by an estimation of the amount of entangled 
mode pairs in momentum space Fig. 2 shows that what 
really limit ates the entanglement is £k,-k-\-q- The bigger 
the entanglement the closest £k,-k-\-q to must be. It 
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FIG. 3. Contour plot of Sk,-k+q for 2J/A = 1.2 and a lattice 
of A/" = 10 sites. Black zones stand for 8 >1 . 



turns out that £k -k+q = iff 2 J/ A > 1 [Cf. below Eq. 
(8)]. 

Criticality and phase transitions. — We have found 
that entanglement diverges in the limit rj ^ 1 when 
£k-k+q = 0. In this limit it is easy to show that p* 
becomes a product of two-mode EPR state [36] , with di- 
verging correlation length. This is an example of a critical 
point which is entirely driven by dissipation [16]. 

However, criticality in dissipative systems is character- 
ized by the vanishing of the /^-eigenvalue with the largest 
real part (closest to zero), Re(Ai) [37]. A practical 
consequence is that the relaxation time becomes infinite, 
Ti ~ [Re(Ai)]~^ oo, close to the transition point. In 
our case Ti ~ (1 — 7^)~^. Consequently the rise in the 
time needed to reach the stationary entangled state needs 
to be considered in any practical application. In other 
words, there is a tradeoff between the maximum entan- 
glement achieved and the time to get it, as exemplified 
by the numerical simulations in Fig. 4. 

Quantum light emission. — So far, we have consid- 
ered the two mode entanglement in the coupled cavity 
system. The question now is how to extract this entangle- 
ment. In doing so, we study the architecture depicted in 
Fig. lb. The resonators are surrounded by an open trans- 
mission line. The line-resonators coupling can be written 
diS H = i Y^7/27r / {e'^^^ ajb{uj) — h.c), where the new 
parameter 7 accounts for the line-resonators coupling 
strength per mode. In the limit of weak coupling, the 
emitted light may be calculated using the input-output 
formalism [33, 38]. We thus obtain 

6rw = ^rw+2y^e-^w, (11) 

where 6|^(t) = l/27r jj^^^ duje-'^%{uj) and 6 = ir/N. 
This is the final result of the letter. By considering the 
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FIG. 4. Time evolution of the Logarithmic Negativity for 
8k,-k+q = and different values 77. 



extended coupling the resonator field modes are mapped 
onto the line with a resolution in frequency 1/A^ [N the 
number of oscillators). Therefore the system presented 
here works as a source of chiral quantum light. 

Conclusions. — In this work we present a practi- 
cal application of a recent paradigm on non-equilibrium 
many body physics: dissipation driven criticality. More 
precisely, we show how to use dissipation and collective 
phenomena as a resource for a source of chiral two-mode 
entangled light. 

From a practical side, our proposal uses simple ele- 
ments, such as qubits with a tunable gap [39, 40], cou- 
pled resonators arrays [23, 24] and external drivings, 
which are standard in state-of-the-art circuit QED ex- 
periments. We argue that it could be thus implemented 
and tested using recent advances in the field of microwave 
state tomography [11, 41, 42]. Our entangled light source 
also features interesting properties: it is directional, with 
entangled modes that can be spatially separated; it is 
broadband, entangling different modes and frequencies, 
and the only parameters that need to be tuned are the 
intensity and the frequency of the external drivings. 

From a theoretical side, our proposal shows that cir- 
cuit QED is a feasible architecture for observing dissi- 
pative many-body phenomena, with concrete examples 
on how to implement these ideas in the lab. We believe 
that the models presented here can be easily extended to 
include novel ingredients, such as photon-photon inter- 
actions and blockade mechanisms, which would open the 
door to an even large family of models, now outside the 
real of Gaussian states. 
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under projects FIS2009-10061 and FIS201 1-25167, by the 
Aragon (Grupo FENOL), QUITEMAD S2009-ESP-1594, 
and the EU project PROMISCE. DP is supported with 
RyC Contract No. Y200200074. 
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APPENDIX A: ADIABATIC ELIMINATION AND 
THE MASTER EQUATION 

In this appendix we detail how to obtain the master 
equation (6). We start with manipulating the coherent 
part. We rewrite Hqq in Eq. (1): 

HcQ = Ho ^ He ^ Hdrive (Al) 

with 

^0 = (i^i + ^ra]aj^ (A2) 
j 

He = XI (4 + %) + ^(4 + + ^i+i) 

(A3) 

Hdrive =^^^cx COs{nat + (A4) 

j a 

where the site-dependent driving phases are given by 

Expanding the resonator operators in momentum 
space (plane wave basis), ak = N~^^'^'^-e~^^^aj^ with 
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k G 27r/N X Z, we can rewrite the total Hamiltonian as 

HcQ = Hq ^ H'^ ^ Hdrive (A5) 

where 

^o = Yl + ^k4ak (A6) 

j k 

^c = E^?(e"''M4 + h.c.) (A7) 

and Sk = uJr-\-2J cos{k), which is vahd whenever J ^oor- 
In the interaction picture with respect to Hq^H driven the 
interaction Hamiltonian is written as 
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{e'^' (a+e2*/(*) + aje-^'^^'^) afee"*^^* + h.c.} 

(A8) 



where the time-dependent term is given by 

f(t) = I* + E E ^ + '^"^>i (A9) 

j a ^ 

as can easily be obtained by integration. We will select 
two driving frequencies = e — {ur + A) and 1^2 = 
e + {ujr + A). Finally, neglecting those terms rotating 
with ujr we get 



H J {t) = Y{a+cj{t) + h.c.) 



(AlO) 



with Cj{t) given by 



+^26" 



-i(/c+02bpi(2Jcos/c-A)t t 



<) (All) 



The couplings terms gi and ^2 depend on the frequen- 
cies 1^1 and 1^2 as 



2Ai ^ 


) Jo 


e - f^ly 


2Ai ^ 




e - l^ly 



2Ao 



: + ^^2 
2A2 



^1 = gJi 
92 = ^-^0 

We will now proceed with the master equation 
^ = Lo(e) + Li(e). 



(A12) 
(A13) 

(A14) 



Here Lq describes the dissipation on the qubits induced 
by the bath - recall that we only take into account spon- 
taneous emission processes, therefore 

p 



while Li describes the Hamiltonian evolution of the cou- 
pling {LiQ=-i[H'^{t),Q\). 



As we want to study the dissipative dynamics induced 
on the resonators by the qubits with steady state |0)(0|, 
we can eliminate adiabatically the irrelevant degrees of 
freedom. We start by defining the projector P 

Pq = /J^^Qq^ss = fJ^^ Qql,ss^ '-'^ Qqi,ss'-'^ QqN,ss' (Al6) 

Here /i describes the system of resonators and we take the 
steady state of the qubits Qqi^ss = as a fixed state 

for them. In perturbation theory up to second order in 
Li we get that 



Itt 



Jo 

(A17) 

where the Born-Markov approximation has already been 
performed. Expanding the commutators in A 17 and tak- 
ing into account that and aj are eigenstates of the 
super operator Lq both with eigenvalue — r/2 will show 
that the only dependence on the integration variable r 
comes through terms of the form exp(— (r/2 ± iujk)r), 
where uj^ = 2 J cos k — A. Integrating these terms we get 



2/r 



i + 2icjfe/r 



(A18) 



According to the hierarchy of energies considered in 
this work, 

e > r > cj^,^, J, A , 

we will always have that Uk/T <C 1 and therefore we 
can neglect the imaginary term in the denominator. Ex- 
panding the commutators in A 17 also shows that the 
associated time evolution of the operators is given 
by a/e exp(— z(2Jcos /c — A)t). Going backwards to the 
Schrodinger picture implies cancelling out these rotating 
terms. We can fulfill this condition applying the following 
transformation 



ak ^ Uo{t)aU^o{t) 



(A19) 



with Uo{t) = e-^^o^ to both sides of A17. Here 0{t) 
denotes an operator in the interaction picture. Summing 
over the sites j we finally get the quantum master equa- 
tion for our coupled set of resonators 



^ = ^-icjfc[4afc,/i] + Y^2bkgbl - {blbk,^}) (A20) 



where bk = giak + 92a^-k+q' 
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Appendix B: Covariance matrix and logarithmic 
negativity 

The QME in momentum space (6) yields the time evo- 
lution for the correlators: 



d_ 



{akai) = -i{uk + uji) - -{gl - gl) {akai) 



-.gig2^i-k+q (Bi) 



dt 



+fg2hi 



(B2) 



and the conjugate equation for (alaj). Both the time 
evolution and the stationary state can be obtained. Note 
that if the momenta k and / do not satisfy the relation 
/ = —k-\-q^ the covariance matrix will be diagonal and 
therefore there will be no entanglement between the emit- 
ted photons. 

If £k,-k+q = 0, the time-dependent covariance matrix 
will take the simple form 



(B3) 



fa -b 0\ 

a b 

-b a 

\0 b aj 

where a and b are time-dependent functions defined by 



a{t) 
b(t) 



gt 



-gl 



^{gl-gl) 
gig2 



m -1+1 



gf 



(B4) 
(B5) 



with F{t) = 1 - exp (- 
negativity is defined as 



r(^i - ^i)^)- The logarithmic 



I 4 

^iV = -^^log2[min(l,2|/,| 



(B6) 



where li are the symplectic eigenvalues of the covari- 
ance matrix. The symplectic spectrum corresponds to 
the eigenvalues of the matrix ^^7^^;.^^, where ft is the 
symplectic matrix 



/ 1 0\ 

-10 

1 
\ 0-10/ 



(B7) 



and the superindex T4 denotes the partial transpose. 
This can be obtained by the transformation 7fc^^/c+g ~ 
A^jk,-k-\-q^ with A defined as 



A 



/I \ 
10 
10 

Vo -1/ 



(B8) 



From definition B6 and matrix B3 we finally get 



EN{t) = -log2 



1+7? 



F{t) + 1 - F{t) 



(B9) 



{v = g2/gi)- 

Following a similar approach for the general case 
£k,-k-\-q 7^ we get the steady state relation 10. 

Appendix C: Input-output formalism 

We consider here the interaction of the coupled qubit- 
cavity (1) and an infinite transmission line (TL). The 
total Hamiltonian for this system is: 



H = Hi 



CQ 



TL 



(CI) 



here i^cQ stands for (1) and i^xL describes the transmis- 
sion line 



TL 



/oo 
dojub^ ( 
-00 



uj)b{uj) 



(C2) 



(where the b operators satisfy the commutation relation 
[b{uj'),b\uj)] = 8{uj' — uj)) and i^int is the interaction 
Hamiltonian - expanding the aj operators in momentum 
space (plane wave basis) it is given by 



jk 



duj 



(^e'^^-k^3 alb{uj) - h.c.) 



Approximating the sum 



(C3) 



to a rectangle of height N centered at k with a width of 
25 — 2i\ jN (being zero elsewhere), we can rewrite C3 as 



pk+6 

V{ / dujNg{uj)alb{Lj) - h.c.} (C4) 

u Jk-5 



From CI, C2 and C4 the Heisenberg equation of mo- 
tion for the operator b{uj) reads. 



b{uj) = —iujbiuj) — \/ ^N/27iak . 



(C5) 



Here, the first Markov approximation giuj) = g = \J^l2i\ 
was assumed. Integrating C5 from to t {to < t) gives 

b{u) = e-'^^'-''^bo{uj) - ^/W/^ f dt'e-'^^'-''^ak{t') 

J to 

(C6) 

Integrating now C6 over the frequency interval [k — d^k-\- 
6] we obtain 

dujb{uj) = V2^bin{t) - [ dt' 



to 



^^-iKt-t') _ t'))ak{t') (C7) 
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with bin defined as (following [38]) 



27r 



k-6 



duje-''^^^-''^bo{uj) (C8) 



In a similar fashion we can integrate C5 from t to ti 
{t < ti) and define a corresponding out operator. We 
will find that the in and out operators are related by 



In the continuum limit ^ 0), taking to 
for sufficiently long values of t, C9 yields 



bin{t) =bout{t) ^ ^1 '^^ak 



with a/c defined by 



-oo and 



(CIO) 



bin{t) = bout{t) + 



dt'- 



y^^-ikit-t') _ t'))ak{t') (C9) 



flfc = ^ / dt'e'^^' ak{t') 



(Cll) 



